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Surfaces of index of refraction have been calculated for the four modes which exist in a warm 
plasma, ignoring all damping mechanisms. The surfaces are displayed as contour diagrams. Dis- 
persion curves are obtained as profile cuts through the surfaces. 



1. Introduction 

The linearized two-fluid plasma theory gives four kinds of waves (four modes); i.e., for every 
combination of plasma parameters there are four independent plane wave solutions to the linear- 
ized set of Maxwell and Euler equations. Many of the properties of these modes have been 
summarized in various books [e.g., Denisse and Delcroix, 1961; Stix, 1962; and Allis, Buchsbaum, 
and Bers, 1963]. The new material in this paper consists of details of the surfaces of n, the index 
of refraction, along with techniques for separating the four modes. 

The plasma parameters are (subscripts e and i are added later when appropriate) 

m = mass of a particle, 
q = charge on a particle, 
N = particle number density, 
T= plasma temperature, 
B = static magnetic field strength, 
= propagation angle, 
co = wave frequency. 

We shall, however, generally use the following normalized parameters: 



* = (a>2 + ft2 )/a> 2 



u = n 2 

a = c 2 /V? 
P = c 2 /V 2 
y = c 2 \V 2 s 
K=cos 2 
S = sin 2 



1 Portions of this work were submitted by Tyan Yeh to Cornell University in partial fulfillment of the requirements for the degree M.Sc. in September, 1963. 

2 Department of Applied Mathematics, California Institute of Technology, Pasadena, Calif. 
:) Center for Radiophysics and Space Research, Cornell University, Ithaca, N.Y. 
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where 

u>p=(Neq 2 Jeome) 112 
n p = (yVigf/e m i ) 1/2 
ft)6=— q e Blm e 
£l b = qiBjmi 
V\ = iKTIme 
V} = 3KTIm 

c = velocity of light 
e = permittivity of vacuum 
k = Boltzmann's constant. 

We shall regard this system of parameters as consisting of the three variables x, y, and u; and 
the four independent parameters Z, r, a, K where Z = qil(—q e ), r=miKZm e ). Note that cojj/ftjj 
= mIClb = r,Vt/Vf = fi/a = rZ, y/a = (1 + r)/(l + 1/Z) and also that 1< a < y < j8. 

For convenience in comparing this work with that in the books mentioned above, we display 
some equivalent symbols in the following table: 



Cohen and 


Denisse and 


Stix 


Allis a et al. 


Yeh 


Delcroix 






u 


y 


ri 2 


ri 1 


X 


X 


a 


a 2 


y 


x/A 


oily 


/3 2 = 0_/3 + 



a Minor modifications are introduced when Z/l. 

The x-y plane is extensively used by Allis et al. [1963] as a background for the CMA diagram 
(so-named by Stix [1962]), which displays phase-velocity surfaces in 13 characteristic regions of the 
x-y plane. These surfaces show velocity as a function of 6, and are to be distinguished from the 
surfaces shown later in this paper, which are contour diagrams of n 2 for a constant 6. Simplified 
versions of these contour diagrams have been used in radio astronomy for many years. [Pawsey 
and Bracewell, 1955.] 

2. Dispersion Equation 

The linearized two-fluid theory of warm plasma gives the dispersion equation as the vanishing 
of a ninth order determinant. Landau damping does not appear in the theory, even though it is 
dominant in some regimes. Collision damping may be included, although we do not do so here. 
The determinant, as given by Denisse and Delcroix [1961], may be expanded to give the dispersion 
equation in the following form: 



Ay 2 -By+C = 
540 



(1) 



with 

A=A\x -\-Aq 

B = B 2 x 2 + B l x + B< ) 

C = C^x* + C2X 2 -h C\x -f Co 
where A = (a - l)(Ka - l)(^/y - 1) 

A = (u - l) 2 (Ka/a - l)(/fo/j8 - 1) 

B 2 = (u- l)(Kuly - 1) + (Kii - l)(u/y - 1) 

Si = (tt - l){(u/a - lX&/i8 - 1) + (Ka/a - l)(a/j3 - 1) 

+ (r+Ur)(u-l)(Kuly-l) 

+ Su[(r+l/r) 

+ a-Hl-l/Z) + )8-Hl-Z)-(l + l/y)]} 
Bo = (a- l) 2 [l/r(a/a - l)(/fo/j3 - 1) + rtXtt/a - l)(a//3 - 1)] 
C 3 = uly — 1 
C 2 -(a/a-l)(u//3-l) + 2(a-l)(a/y-l) 

Ci = (u-l)[2(u/a-lXtt/j8-l) + (a-lXa/y-l)] 

Co = (u-mula-l)(ul(3-l); 
hence ^ =( u - l)[(K tt - l)(Kuly- l)x + (u- l)(Ku/a- l)(Kulf3- 1)] 

B = K(u-l)[2x + (fila + aiP)(u-l)] [(u/y-l)x 

+ (a/o-lXtt/j8-l)] 

-S(* + a - l)[(a + a/7 - 2)x 
+ (a - lXa/a - 0/a - n/j8 - a/0)] 
C = (x + a-l) 2 [(a/y-l)^ + (Wa:-l)(a/i3-l)]. 

3. Some Features of the Contours 

Equation (1) is a quartic in u. Denisse and Delcroix [1961] have shown that all four values 
of u are real. In the x-y-u space, (1) is a surface of refractive index, consisting of four sheets, 
one for each mode. For a given value of u, (1) becomes a contour on the x-y plane. We are 
interested in the structure of the surface of refractive index, but only in the first quadrant. 

Some features of the surface can be found readily in analytic form. These are: 

(1) The behavior at the origin. 

(2) The contour u = 0. 

(3) The contour u=l. 
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(4) The contours u = ± c °. 

(5) A contour possesses two x-axis intercepts (at one, the contour is tangent to the x-axis; at 

the other, it crosses the x-axis). 

(6) A contour possesses two y-axis intercepts (it crosses the y-axis at these two intercepts). 

(7) A contour possesses two oblique asymptotes. 

(8) A contour possesses one vertical asymptote. 

(1) At the origin: When x = and y=0, (1) becomes 

(u-iy(u/a-l)(ulp-l) 

the four values of u are 1, 1, a, and p. 

(2) The contour for u = (fig. 1): When u = Q, (1) becomes 

(l-*)[(y + *-l) 2 -&y] = 0, 

where b = (r 1/2 — r~ 1/2 ) 2 . This contour consists of the line x=l and a parabola whose axis is in- 
clined to the .%-axis by 37r/4 and whose vertex is (1 + 36/16, 6/16). The parabola is tangent to the 
x-axis at (1, 0). 

(3) The contour for u— 1 (fig. 2): When u= 1, (1) degenerates to 

x 2 [S(l/y-l)y + (l/y-l)x + (l/a-l)(l//3-l)] = 0. 

The contour consists of two vertical lines % = 0, and one oblique line with slope —1/S. 

(4) The contour for u = ± cc (figs. 3, 4): When w=± 00 , the coefficient of u 4 in (1) vanishes, 
giving 

(llap)(Ky-r)(Ky-llr) = 0. 




HI 



/3- 



/3-K 



, ( i- ,)( /r'> 



I a-l 



r a-K 



m 



m \ 



<>' 



\nz: 

\ 

\ 



ii 



Figure 1. The contour u = 0. 
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Figure 2. The contour u = l. 
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Figure 3. 77ie contour u = + <*>. 



Figure 4. 77ie contour u = — ». 
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FIGURE 5. Intercepts of the contours on the x-axis, as 
functions of u. 



The contour consists of two horizontal lines y = r/K and y=l/Kr. In addition, the contour for 
u= — °° has a vertical line at x= + °° (see (7) below). 

(5) The x-axis intercepts (fig. 5): When y=0, (1) becomes 

fee + tt - l) 2 [(u/y - l)x + (ula - 1) (u/jS - 1)] =0. 

At the double zero % = ^(a) where 

*,(a) = 1 — u (2) 
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the contour is tangent to the x-axis. As the simple zero x = x c (u) where 

x c (u) — (u/o- 1) (tt//3-l)/(u/y-l) 



(3) 



the contour crosses the %-axis. 

(6) The y-axis intercepts (fig. 6): When x = 0, (1) becomes (u—\) 2 [{Kuja— l)y—r~ 1 (ul<x— 1)] 
[(Ku/13— l)y— r(«//3 — 1)] = 0. At the two zeros y = y a (u) and y = yb{u) where 



y a (u) = 7W/3-l)/(JWj3-l) 



(4) 
(5) 



the contour crosses the y-axis. The double zero u= 1 is part of the contour described in (3) above. 
(7) The oblique asymptotes (fig. 7): When x — > °° while yjx is fixed, (1) approaches 

x[(u- 1) (Ka/y- l)y-(u/7~ 1)%] [(Ka- l)y-x] =0 
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Figure 6. Intercepts of the contours on the y-axis, as 
functions of u. 
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FIGURE 7. Slopes of the oblique asymptotes to the con- 
tours, as functions ofu. 
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FIGURE 8. Positions of the vertical asymptotes to the con- 
tours, as functions of u. 



which gives two oblique asymptotes y = S\(u) x and y = Sz(u)x where 

si(tt) - (uly- l)l(u-l)(Ku/y- 1) 
s 2 (u) = ll(Ku-l). 



(6) 

(7) 



Now (1) is a quartic in u, and the product of the four roots is — a/3(x— l)[(y + x— l) 2 — (r 1/2 — r _1/2 ) 2 y]/ 
(Ky— r)(Ky— 1/r, the ratio of the coefficients of u° and u 4 in (1). As x — ► », one root is (1 + xly)/K, 
given by (7), and the product of two other roots is y(l-\- x/y)/K, given by (6). Hence the fourth 
root is — a(3(x — l)/y, which shows that there is a contour curve running vertically at x = + <x> when 

(8) The vertical asymptote (fig. 8): When A vanishes, which occurs at x=x O0 (u), where 



xM ^-(u- l)(Kula - l)(Ku/p - l)l(Ku - l)(Ku/y - 1) 



(8) 



either the root y+ (if B > at x = x cc ) or the other root r_ (if B > at x = x O0 ) becomes infinite [see 
(9) and (10)]. In the vicinity of the vertical asymptote the behavior of the contour is determined 
by the sign of A\. 



4. Contour Curves 

By the contour of a given u, we refer to the two branches in the first quadrant of the x-y plane 
given by y = y+(x) and y = y~(x), 

where 



and 



y+(x, u) i 



y-(x, u) 



b±Vd 

1A 

b-Vd 

2A 



(9) 
(10) 
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where D = B 2 — 4>AC. By a contour curve we refer to one continuous section of a contour; it starts 
at a point on one of the axes, or on an asymptote, and terminates on the other axis, or along an 
asymptote. Depending on u, a contour curve can be given solely by the branch y = y+{x, u), by 
the branch y = y~(x 9 u), or by a combination of the two branches joined smoothly at an extreme 
point x = x e (u) where D vanishes. There can be one, two, three, or four contour curves in a contour. 

When D < 0, the contour does not exist. When D = 0, the two branches join smoothly at 
the extreme point. When D > 0, the two branches are separated. The two branches never cross 
each other, at most they can be tangent to each other at an extreme point. For in the vicinity 
of an extreme point, D goes through zero (a) from negative to negative, or (b) from positive to 
positive, or (c) from one sign to the other. In case (a) the extreme point is an isolated point. In 
case (b) y+(x) ^ y-(x) if A > at x = x e or y+(x) ^ y-(x) if A < at x = x e . The case A = at x = x e 
is trivial because then all A, B, and C vanish. In case (c) the contour exists only to one side of 
the extreme point. 

The fact that the two branches y+(x, u) and y~(x, u) can not cross each other implies that the 
four sheets of the refractive index surface do not penetrate into each other. For, each sheet comes 
either from y+(x, u) or from y~(x, u); if they interpenetrated at some u = Uq the corresponding contour 
curves for u = Uo would cross, but they do not. 

In the remainder of this section we discuss in detail the special features of section 3, and iden- 
tify the contour curves with y+ and y_. In section 4 this identification will be used to help sort 
the curves into the four modes. 

In the figures which follow, the solid curves result from the branch y = y+(x), and the dotted 
curves result from y = y~(x). 

(1) The contour for u = (fig. 1): When u = 0, 

A=-{x-\) 

B = (x-l)[2x-(r+l/r)] 

C=-{x-lf 

D=-b{x-lf{^x-a) 

where a = (r 1/2 + r~ 1/2 ) 2 and b = (r 1/2 — r~ 1/2 ) 2 . The parabola may be split into two parts: the upper 



half yi(x)=-x + i~\a + b) + 2~W b(a - 4*) and the lower half y 2 (x)=-x + ^~\a + b) - 2" 1 Vb(ar-4x). 

The two halves join at the extreme point (a/4, 6/4). Since VD = — (x — 1) V6(a — 4>x) for x ^ 1 and 

VD = (x— 1) V6(a — 4*) for 1 ^ * ^ a/4, then y+, - = yi, 2 (x)forx ^ l,andy+, -=y 2 , i(x) for l^x^a/4. 
The appropriate sections of the parabola in figure 1 are solid and dashed, corresponding to y+ and 
y_. The vertical line x= 1 is split up as shown in figure 1 by consideration of the limit u^> 0. 

(2) The contour for &= 1 (fig. 2): As w^l, we have, to first order, 

A^-S(u-l)(K/y-l)x 
B->-S(l/y-l)x 2 

C^x 2 [(l/y-l)x + (l/a-l)(l//3-l)] 

Z)->S 2 (l/y-l) 2 ^ 4 -f4S(a-lXX/y-l)x 3 [(l/y-l)^-h(l/a-l)(l/i8-l)] 
VD-^-S(lly-l)x 2 -2(u-l)(K/y-l)x[x^(lla-l)(ll/3-)l(lly-l)]. 

Hence y_^-S- 1 [^ + (l/a- l)(l/j8- l)/(l/y- 1)] and y+-> (1/y- l)x/(u- l)(K/y- 1), and the 
oblique line is given by y~(x\ as shown in figure 2. A consideration of the limit u — > 1 shows that 
one of the vertical fines is given by y+, and the other vertical line is given by y_. 
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(3) The contour for u = ±°° (figs. 3, 4): Asm~^±oo, 

A-*K 2 u 4 laP 

B^(r+l/r)Ku 4 laP 

C^>u 4 lap 

D^(r-l/r) 2 K 2 u 8 /a 2 p 2 . 

Since V5-*(r- l/r)Ku 4 /ap, then y + -* r/K and y. -» UKr. 

(4) The *-axis intercepts (fig. 5): When x = x t (u) 9 B = -Kbu{u- l) 2 [(l/y- l/a/8)u-(l + 1/y) 
+ (1/a + l/j8)l which is non-negative if and only if < u < (1 + 1/y — 1/a — l//3)/(l/y — l/a/3). Hence 
the tangency point results from y+ when a < 0, and from y_ when < a < 1. When % = x c (a), 
B=-Su[(l/-y- l/a/8)tt~(l + 1/y- l/a-l/j8)]/a)/(a/y- l) 2 where /(a) - (a- l)(a/y- DKr^X^/a 
— 1) + r(^//3 — 1)] — (&/a — l)(u//3 — \){u — 1 + u/y — 1)] , which has only 4 one real zero, denoted by 
u = 8, lying between and 1. Since B is negative if and only if < a < 8 or u > (1 4- 1/y— 1/a 
— l//3)/(l/y— l/a/3), the crossing point results from y+(u) when < a < 8 as well as y < a < ft and 
from y _ [u) when & < as well as 8 < u < a. 

(5) The y-axis intercepts (fig. 6): When x = 0, 

A = (u- lf(Ku/a - 1) (Ku/p - 1) 

B = (u - l) 2 [r(Ku/a - 1) (a/jB - 1) + r-^w/a - 1) (Ku//3 - 1)] 

C = (a-l) 2 (tt/a-l) (*//£- 1) 

Z) = (u-l) 4 [g(")] 2 

where g(u) = r(Kula — l) (u//3—l) — r- l (ula—l) (Ku/f3—l). The discriminator of g(u) may be 
written 

ifilZ -aZ)- 2 {(r-l/r) 2 -[VrZ - 1/ V72) 2 } 2 {K- K X ){K- K 2 \ 
which is negative if and only if K is in the very small range K\<K< K 2 , where 

Ki = ( V7/Z-VzTr) 2 l(r- 1/r + V7Z - 1/ VTZ) 2 

K 2 = ( V7JZ- VzTr) 2 l(r- 1/r- Vrf + 1/ VTZ) 2 . 

If this inequality holds, g(u) is positive-definite and the two curves y a (u) and yb(u) [in (4) and (5)] do 
not intersect. Since vD = (u— l)?g(u), y+ = y«(w) and y_ = yb(u). Otherwise, g(u) is negative if and 
only if u\ < u < u 2 where u\ and u 2 are the two real zeros of g(u), lying between a/K and j8. The two 
curves y a {u) and yb{u) intersect at u\ and u 2 , and y+ and y_ switch between y a (u) and y&(a) at &i and 
u 2 . Thus, for u\<u< u 2 , vD = — (u — l) 2 g(u). Hence y+ = yb{u) and y-=y a (u). For u < u\ as 
well as u > 1/2, \D=(u— l) 2 g(u), hence y+ = y a (u) and y_ = y&(^). 

Figure 6 is drawn for the case K 2 < K^ 1; i.e., not too close to 90°. The curve y = y a (u) 
and the line a = \ intersect at the point y=ri(3— 1) (ft — K). The curve y = yb{u) and the line 
u=\ intersect at the point y=r~ l {a— 1) (a — K) . 



4 This point has been established by numerical calculations. 
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When K = K 2 , where a/K < u\ = u 2 < /3, the two curves are tangent at a point in the first 
quadrant. When K X <K<K 2 the two curves are separate. When K = K U where P<ui = u 2 
< a/K, the curves are tangent in the fourth quadrant. When < K < Ki the curves intersect 
in the fourth quadrant. Since only positive y has physical significance, we can state that the 
curves switch between y+ and y_ unless K < K 2 , or, approximately, 9 > cos -1 (m e lmi) 112 . 

(6) The oblique asymptotes (fig. 7): As *— »°° while y\x is fixed, 

A-+(u-l) (Ku-l) (Kuly-l)x 
B-+[(u-l)(Kuly-l) + (Ku-l)(u/y-l)})c 2 
C^(u/y-l)x* 
D->[h(u)] 2 x 4 

where h(u) = (Ku- 1) (u/y- l)-(u- 1) (Ku/y- l) = Su(l - 1/y). For u 2* 0, h(u) is nonnega- 
tive, and V r D = h(u), y+(x, u) = s 1 (u)x, and y~(x, u) = s 2 (u)x. 

(7) At the origin: When u = a and x = 0, A = -S(a~ l) 2 (Kalp- 1), B = S(a- 1) 2 (1- a//3)r, 
C = 0,V5 = S(a-l) 2 (l-/3/a)r; hence y_ = 0. When n = j8 andx = (M = -S(j8- l) 2 (Kp/a- 1), 
5 = -S( i 8-l) 2 (i8/a-l)/r,C = 0,V5-S( i Q-l) 2 (/3/a-l)/r; hence y + = 0. 

(8) A further special case which will be used later is the calculation of y+ when u = y and 
x^+°o: In this case A^-S(y-l) (Ky-l)x, B-+-S(y-l)x 2 , C-» (y/a- 1) (y/j8- l)^ 2 , 
V5 -+-B + 2AC/B; hence y+^C/fi=- (y/a- 1) (y/j8-l)/S(y- 1). 



5. Separation of Modes 

The foregoing analyses allow us now to separate the four modes; that is, to identify contour 
curves for various values of u with one of the four sheets of the index of refraction surface. This 
separation can be done in an unambiguous way because the four sheets do not penetrate into each 
other, and the contour curves do not cross, as shown in section 4. 

Consider figures 5, 6, and 7. Each is a profile of the index surface, made by cutting the sur- 
face by the planes y = 0, x = 0, and x=-\-°c, respectively. The curves labeled I, II, III, IV must 
each belong to one mode, else (in figs. 5 and 6) the curves would cross. We conjecture that the 
introduction of collisions would remove the points of contact and would leave the curves entirely 
distinct as they already are in figure 7. 

Curves 5III and 6III have a common point (x y y, u) = (0, 0, /3). They both come from y+(#, u) 
and hence are different profiles of the same sheet of the index surface; they belong to the same 
mode. Curve 5III is symptotic to the point [+ 00 , — (y/a— 1) (y/fi— 1)/S(y— 1), y], which cor- 
responds to the point (0, y) on curve 7III. Curves 5III and 7III both come from y+(x, w), and 
these belong to the same mode. The three curves 5III, 6III, and 7III therefore all belong to one 
mode, call it mode III. By the principle of nonpenetration, the three curves 51 V, 61 V, and 71 V, 
immediately to the left of 5III, 6III, and 7III, must be profiles of one mode; call it mode IV. Simi- 
larly, the curves labeled I and II must belong to the other two modes. 

Any contour curve at least has a point on the x-axis or the y-axis, or is asymptotic to an oblique 
asymptote. From the identification of these features with modes according to figures 5, 6, and 7, 
the mode to which a contour belongs is uniquely determined. Ambiguities which arise at one end 
of a contour (e.g., the contour which has a y-axis intercept at the point of contact between curves 
611 and 6III) are resolved at the other end of the contour. 

In this fashion we can construct the four index surfaces in a unique way. If we look from 
above towards the origin, we see (in the first octant) the sheets in the sequence I, II, III, IV, then 
two further sheets of modes I and II. This nomenclature agrees with that of Denisse and Delcroix 
[1961]. This method of separation is valid for any propagation angle. The details of the surface 
change with 6, however, especially as — > 90°. 
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6. Index of Refraction Surfaces 

Figures 9, 10, 11, and 12 show, in a schematic fashion, the index surfaces for the four modes; 
the solid and dashed portions of the curves correspond to y+(x, u) and y~(x, a), respectively. 

Figures 13, 14, 15, and 16 show the detailed surfaces for the case r= 100, Z = l, a = 50, = 30°. 
This is not a realistic plasma, but this choice of parameters does display the various features of 
the surfaces in a convenient manner. The contour curves are labeled with their values of u. 
The calculations are made from (1) by choosing sets of values of u and x and solving for y on a 
digital computer. The resulting curves are sorted into modes by examining endpoints and asymp- 
totes, as described in the preceding section. 
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Figure 9. Schematic contour diagram of index of refrac- 
tion for mode I, on x-y plane. 



FIGURE 11. Schematic contour diagram of index of refrac- 
tion for mode III, on x-y plane. 
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FIGURE 10. Schematic contour diagram of index of refrac- 
tion for mode II, on x-y plane. 



FIGURE 12. Schematic contour diagram of index of refrac- 
tion for mode IV, on x-y plane. 
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Figure 13. Contour diagram of index of refraction for mode I, for r — 100, Z — 1 , a — 50, 6 — 30° 



Curves are labeled with numerical value of u = n 2 . 




FIGURE 14. Contour diagram of index of refraction for mode II, for r = 100, Z = 1, ot = 50, 6 = 30° 



Curves are labeled with numerical value of u = n 2 . 



550 




\(f)y(\ + \0'°)y 



FIGURE 15. Contour diagram of index of refraction for mode III, for r = 100, Z = 

Curves are labeled with numerical value of u = n 2 . 
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FIGURE 16. Contour diagram of index of refraction for mode IV, for r = 100, Z = 1, a = 50, 6 = 30° 

Curves are labeled with numerical value of u — n 2 . 
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We observe that 

(1) Mode III has u > for all x, y, whereas modes I, II, IV have u < in some regions of the 
x-y plane. These are regions where the corresponding waves are evanescent. 

(2) Modes I, II, IV contain the line w = 0, which is the boundary between propagation and 
evanescence. A point u = is called a cutoff. 

(3) Modes I and II contain singularities, u = ± (X > . A point w = + °°is called a resonance. 

(4) Modes II, III, IV contain regions where the surface has a very steep slope, falling abruptly 
from yjK to a/K (mode II), from a/K to unity and from /3 to a (mode III), and from a to unity (mode 
IV). These regions are called pseudo resonances. In the limit 7 1 — » they turn into the true 
resonances of magnetoionic theory. [Denisse and Delcroix, 1961.] 

(5) The regions of oblique asymptotes of modes I, II, and III are the regions of the three mag- 
netohydrodynamic waves. Modes I and III give the slow and fast magnetoacoustic waves, while 
mode II gives the Alfven wave. Equation (7) gives the Alfven velocity, divided by K, as the well- 
known asymptotic phase velocity of mode II. 

(6) The valleys in modes II and III give minima in the dispersion curves. The valley of mode 
III is the region of common whistlers. The valley of mode II is a region of MHD whistlers. These 
MHD whistlers do not appear to have been discussed much in the literature, although the minimum 
in the mode II dispersion curve shows in figure VI, 10 of Denisse and Delcroix [1961]. As shown 
in figure 18 below, this minimum does not appear except at very small values of y/x. 

(7) The cold plasma modes can be found by letting a, /3, y — » « . The surfaces can then be 
coalesced into 2 sheets, the ordinary and extraordinary waves. 



7. Dispersion Curves 

Dispersion curves giving u(o)) are obtained by cutting the index of refraction surfaces by planes 
y/x = (ObQb I ((*)% + ft*). Each plane is specified by a specific ratio of plasma density to strength of 
the static magnetic field. The dispersion curves for each mode are confined to limited regions on 
the u-x plane. This is displayed in figures 17 to 20, each of which is the family of dispersion curves 
for one mode, for a wide range of values of (y/x). 

It is of interest to collect together the four dispersion curves for the various modes for constant 
y/x. This is done in figures 21 to 25, which show the combined curves for y/x=l0 4 , 10 2 , 1, 10~ 2 , 
10 4 . 

The dispersion curves are constructed from the index surfaces, figures 13 to 16, and thus apply 
to that particular plasma, although they are similar for other values of the parameters. Denisse 
and Delcroix [1961, See ch. VI] have given a set of dispersion curves comparable to figures 21 to 
25 for the case r=1836, Z = l, a « 10 5 , = 0°. 

In ionospheric work it often is useful to regard the magnetic field as constant while letting the 
density vary. The dispersion curves corresponding to this situation are obtained by cutting 
the index surfaces with planes at constant y. The profiles are again plots of u versus x, but they 
may be interpreted as plots of u versus height. 

In radio astronomy there often is interest in following a monochromatic ray through a region 
where B and N both vary. The set of values (x, y) along the ray's trajectory describes a curve on 
the x-y plane, and the curve can be drawn on each of the four index surfaces. The value and 
gradient of the index can thereby be seen for each of the four waves, and with this information 
one may be able to infer something of the nature of the source. [Pawsey and Bracewell, 1955] 
The simple assumptions concerning the propagation which one generally makes are (a) a wave is 
totally reflected at a cutoff; (b) a wave is totally absorbed at a resonance or at a pseudo-resonance; 
(c) a wave suffers high attenuation if its phase velocity is comparable with the thermal velocity 
of the plasma particles; and (d) energy can be coupled from one wave to another in regions of high 
gradient or in regions where two different modes have nearly the same index. (For theoretical 
discussions of these processes, see, e.g., Stix [1962], Budden [1961], and Ginzburg, [1961].) 
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FIGURE 17. Dispersion curves for mode I, Parametric in y/x = a) b flJ(o) 2 + ft*) , for r= 100, Z= 1, a = 50, 

= 30°. 
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FIGURE 18. Dispersion curves for mode 11, parametric in y/x = ct> b ft b /(aj 2 -|- ft 2 ), for r=100, Z=l, 

a = 50, = 30°. P P 
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Figure 19. Dispersion curves for mode III, parametric in yl\ = (Obflbl(co 2 ^rfl 2 ), for r = 100, Z = l, 

a = 50, = 30°. 
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FIGURE 20. Dispersion curves for mode IV, parametric in y/x = a> b fl b /(a> 2 + n 2 ), for r = 

a = 50, = 30°. 
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Figure 21. Dispersion curves for y/x = co b n b /(a; 2 + H 2 ) = 10 4 ,/or r= 100, Z = 1, a = 50, = 30° 
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FIGURE 22. Dispersion curves for y/x = aj b n b /( w 2 + 2 ) = 10 2 , /or r = 100, Z=l, a = 50, = 30° 
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FIGURE 23. Dispersion curves for y/x = a> b n b /(cuj + nj) = 1, /or r= 100, Z=l, a = 50, = 30°. 
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FIGURE 24. Dispersion curves for y/x = a) b n b /(a)J + nj) = 10" 2 , /or r= 100, Z = l, a = 50, = 30° 
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FIGURE 25. Dispersion curves for y/x = <D b n b /(a>2 + ny = 10 4 , /or r= 100, Z=l, a = 50, = 30°. 
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